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RANDOM GROWTH MODELS WITH POLYGONAL SHAPES^ 
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We consider discrete-time random perturbations of monotone cel- 
lular automata (CA) in two dimensions. Under general conditions, 
we prove the existence of half-space velocities, and then establish the 
validity of the Wulff construction for asymptotic shapes arising from 
finite initial seeds. Such a shape converges to the polygonal invariant 
shape of the corresponding deterministic model as the perturbation 
decreases. In many cases, exact stability is observed. That is, for small 
perturbations, the shapes of the deterministic and random processes 
agree exactly. We give a complete characterization of such cases, and 
show that they are prevalent among threshold growth CA with box 
neighborhood. We also design a nontrivial family of CA in which the 
shape is exactly computable for all values of its probability parame- 
ter. 

1. Introduction. Discrete local models for random growth and deposi- 
tion have been a staple of rigorous research in probability since the Ham- 
mersley and Welsh paper [18] on first passage percolation about 40 years 
ago. Apart from their role as a testing ground for probabilistic techniques, a 
voluminous physics literature [26, 28] testifies to their importance in under- 
standing the evolution of natural systems far from equilibrium. The most 
basic tool, introduced in [18] and ubiquitous ever since, is subadditivity: the 
process dominates one restarted from an already occupied point. Clearly, 
this imposes a monotonicity property on the model, but, as we will see, not 
much more. The result is the existence of an asymptotic shape: started from 
a finite seed, and scaled by time, the occupied set converges to a determin- 
istic convex limit. Elegant as this method is, it is nonconstructive and as a 
result fails to provide any detailed information about the limiting set. Thus 
asymptotic properties of subadditive sequences are still an active area of 
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research [2, 34]. Are there cases when the shape can be exactly identified? 
Research on this topic has so far primarily focused on growth from infinite 
initial states, also known as random interfaces. Methods have ranged from 
hydrodynamic limits based on explicit identification of invariant measures 
[30, 31], to techniques arising from exactly solvable systems in mathemat- 
ical physics [14, 20], and to perturbation arguments based on supercritical 
oriented percolation [6, 8] which imply that some interfaces move with the 
speed of their deterministic counterparts. For other related rigorous and 
empirical results see [16, 22, 27]. 

The main aim of this paper is to extend the perturbation approach to 
show that the finite limit shape of a random growth model may also agree 
with that of a deterministic one. At issue is not merely whether small random 
errors induce small changes (we will see that this is always the case), but 
rather whether the shape can stay exactly the same. This property, which 
we call exact stability, is only valid under substantial assumptions, as the 
model has to have opposite structure, in an appropriate sense, from the 
additive one considered in [6]. In the process, we extend the result of [3] to 
obtain the Wulff characterization of the invariant shape. We also show that 
exact stability is far from rare; in fact, almost all members of arguably the 
most natural family of two-dimensional growth models, the threshold growth 
cellular automata with square neighborhood, are exactly stable. Finally, we 
show how to employ exactly solvable systems to construct one example 
which has a computable shape for every value of its probability parameter. 
(Although they are invaluable in suggesting universal phenomena, exactly 
solvable examples are extremely difficult to come by.) 

The random rules we describe below can be thought of as discrete coun- 
terparts to the KPZ equation, which in turn is touted as a universal scaling 
model for any local growth and deposition process in physics [26, 28], in 
particular for crystal growth [28] . This we mention because the well-studied 
roughening transition in crystallography, whereby a crystal loses its polyg- 
onal shape as the ambient temperature increases, produces pictures which 
are strikingly similar to ours [32]. While this transition is usually thought 
to be an equilibrium phenomenon, the present results at least suggest that 
it may have a dynamic counterpart. 

We now proceed to precise formulations. Unfortunately, these require a 
large number of definitions related to our previous work. Although we do 
not use any of the results from [10] explicitly, a glance at that paper's first 
two sections may help to motivate what follows. 

Our basic framework consists of two-state cellular automata (CA). In 
general, such a cellular automaton is specified by the following two ingre- 
dients. The first is a finite neighborhood M dT? of the origin, its translate 
X + N then being the neighborhood of point x. By convention, we assume 
that J\f contains the origin. Typically, M = By{Q,p) = {x : ||x||^ < p}, where 
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II • Wy is the £'^-norm. When u = \ the resulting N is called the range p Di- 
amond neighborhood, while if i/ = cxo it is referred to as the range p Box 
neighborhood. (In particular, range 1 Diamond and Box neighborhoods are 
also known as von Neumann and Moore neighborhoods, resp.) The second 
ingredient is a map vr : 2^ {0, 1}, which flags the sufficient configurations 
for occupancy. More precisely, for a set A dT?, we define T{A) C 7? by 
adjoining every x G for which '^{{At — x) nJ\f) = 1. Then, for a given ini- 
tial subset Aq C of occupied points, we define Ai,A2,... recursively by 
At+i = T{At). Accordingly, occupied and vacant sites will often be denoted 
by I's and O's, respectively. Our main focus will be starting states which 
consist of a possibly large, but finite set of I's surrounded by O's. However, 
we will also consider other initial states, namely half-spaces and wedges. 

We restrict to two-dimensional dynamics for two main reasons. First, 
almost every step in higher dimensions introduces new technical complica- 
tions, some quite serious. In fact, there are new phenomena, and the classi- 
fication of Theorem 2 below becomes much more complex. Second, some of 
our techniques are intrinsically two-dimensional, such as the explicitly solv- 
able example of Section 6, the lattice geometry and analytic number theory 
of Section 7, and even combinatorial properties studied in [3]. Neverthe- 
less, some results — notably Theorem 1 — do readily generalize to arbitrary 
dimension. 

Our key assumption is that the CA dynamics are monotone (or attrac- 
tive), that is. Si C 52 implies 7r(S'i) < 7r(S'2). Note that specifying a mono- 
tone dynamics is the same as specifying an antichain of subsets of M: the 
inclusion minimal sets S with 7r(5') = 1 having the property that none of 
them is a subset of another. Surprisingly, the number of possible monotone 
dynamics (known as a Dedekind number) is possible to estimate for large M. 
Some typical properties of monotone CA are also known [23] . Unfortunately, 
it turns out that for large box neighborhoods the asymptotic proportion of 
supercritical rules (see the definition below) is negligible. Other interesting 
properties seem to present great difficulties. In studying typical monotone 
CA rules, it is therefore desirable to restrict to a simpler class. 

A natural such class consists of totalistic monotone CA, those for which 
7r(S') depends only on the cardinality 15*1 of S. In other words, there exists a 
threshold 9 >0 such that vr(5') = whenever l^l < ^ and tt{S) = 1 whenever 
15*1 >9. This much studied case is also known by the name threshold growth 
(TG) CA. 

Induced by T is a growth transformation T on closed subsets of , given 

by 

T{B) = {x G G T{{B - x) n Z^)}. 

In words, one translates the lattice so that x G R^ is at the origin, and 
applies T to the intersection of Euclidean set B with the translated lattice. 
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It is easy to verify that the two transformations are conjugate, 

r{Br\7?)=f{B)r\7?. 

It wiU become immediately apparent why T is convenient. Let 5^ C be 
the set of unit vectors and let 

H~ = {x eYi^ : {x,u) <id] 

be the closed half-space with outward normal u G 5^. Then there exists a 
w(u) E R so that 

and consequently 

T*(i7^ n z^) = {H~ + tw{u) ■ u) n z2. 

If w{u) > for every u we call the CA supercritical. A supercritical CA hence 
enlarges every half-space. This is equivalent to existence of a finite set Aq 
which fills space, that is, {J^^q At = T? [3, 10]. All initial sets will be assumed 
to fill space from now on. Set 

Ky^ = [jW,l/w{u)]-u:ueS^} 

and let L be the polar transform of Ki/^, that is, 

L = ^1% = {x£'R^: {x,u) < w{u)}. 
Then one can prove the following limiting shape result for any finite 



where the limit is taken in the Hausdorff metric. In short, the shape L = L['k) 
is obtained as the Wulff transform of the speed function : ^ R, which 
for small neighborhoods is readily computable by hand or by computer. 
Furthermore, L is always a polygon and the Hausdorff distance between At 
and tL is bounded in time t [9, 10, 11, 12, 13, 37]. 

To formulate the stability properties of L under random perturbations, we 
begin by introducing a general monotone random dynamics. The function 
TT differs from the one described above in that it has values in [0, 1] . Upon 
seeing a set of occupied sites a; -|- 5 in its neighborhood at time t, a site 
becomes occupied at time t + 1 independently with probability 7r(5). To 
obtain a monotone rule we require that 7r(5i) < 7r(S'2) whenever C 52. 

More precisely, introduce i.i.d. vectors ^x,t-, x£ T? , t = 0,1,2,..., with 
2l"^l coordinates ^x,t{S), which are Bernoulli(7r(S')) for every S C M. We 
assume that these are coupled so that Cx,tiSi) = 1 implies that ^x,t{S2) = 1 
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whenever Si C 82- The construction of such a couphng is left as an exercise 
for the reader. The random sets Ai,A2,... are now determined by 

At+i = {x : + AA) n At) -x) = 1}. 

To avoid some triviahties and inessential complications, we assume that 
I's only grow by contact: 7r(0) =0, and that vr is symmetric: —N=N and 
7r(— 5) = vr(5). Much more substantial is the assumption that vr solidifies: 
7r{S) = 1 whenever G 5". These three properties, together with monotonic- 
ity, will be our standing assumptions throughout the paper. 

For every random vr, we set p = min{7r(S') : 7r(5) > 0}, define the associated 
deterministic dynamics by its map iTdiS) = l^^(^s)>o}i label the iteration 
transform T as before. We will say that vr is a, p- perturbation of the CA T. For 
many purposes the standard p-perturbation, which has 7r(5) =p whenever 
7r(S') > 0, suffices. 

We say that a p-perturbation of T has shape if 

lim — - = 

t— >oo t 

almost surely, in the Hausdorff metric, for every finite initial set ^0 which fills 
space. We say that T has exactly stable shape L if there exists a p <1 such 
that = L (which of course subsumes the existence of the shape L^r) for 
the standard, and hence any, p-perturbation vr. For a standard perturbation, 
we also write Lp = Lj^. Thus Li = L. Recall that the deterministic growth 
at time t is included in a constant fattening of tLi; hence the same is true 
of any p-perturbation. 

As already mentioned, such considerations are in the general direction of 
the vintage Durrett-Liggett flat edge result [6]. To describe their result in our 
context, recall that a deterministic CA is additive if tt{S) equals 1 precisely 
when S is nonempty. In this case Ki/^ =M* and L = co{M). Moreover, any 
standard perturbation is a first passage percolation model, and as such has 
an almost sure (deterministic) limiting shape Lp for each p > [5, 29]. For 
the von Neumann neighborhood, Durrett and Liggett proved that, if p is 
close to 1, then Lp is close to L and in fact inherits from L flat edges in the 
four diagonal directions. However, they show that Lp is not equal to L, due 
to the fact that its extent in the coordinate directions is strictly less than 1. 

The existence of a limiting shape for general random dynamics does 
not immediately follow from standard subadditivity arguments. A sufficient 
condition is a property of T we call local regularity. Namely, for every initial 
state Aq there exists a constant C so that the following is true for every 
fixed (deterministic) assignment of ^x,t- every x € At at distance at least C 
from ^0 has an occupied set G C At entirely within distance C of x such 
that G fills space. 
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Note that local regularity is a combinatorial condition involving every 
possible way At can evolve, and thus has nothing to do with probability. 
At first it seems a condition not likely to be often satisfied, but the oppo- 
site appears true. One can easily check local regularity directly for many 
cases with small M, and it holds generally for box neighborhood TG CA. 
All known counterexamples involve "strange" neighborhoods [3]. Under this 
condition, it can readily be shown that L^^ exists. 

Besides finite shapes, limiting profiles from half-spaces are of considerable 
interest. The first reason is that their Monte Carlo approximations can be 
computed much more efficiently (see Remark 2 in Section 8). The second is 
that they are important for shapes from other infinite sets, such as wedges 
and holes [13]. For finite seeds also, the Wulff transform (see Corollary 1.1 
below), which expresses the asymptotic shape in terms of half-space veloc- 
ities, is very handy. However, the limit theorem in [3] does not extend to 
infinite seeds, as restarting requires an a priori upper bound on fluctuations. 
Here we provide the missing step, which establishes the following large de- 
viations bound, referred to as the Kesten property in [13]. (See [21] for a 
similar result in the first passage context.) 

Theorem 1. Let tt be a p-perturbation of a locally regular supercritical 
CA and let the initial set be Aq = H~ n for u G . Then there exists a 
deterministic Wt^{u) > such that 

+ t{wT,{u) - e)-u(l AtC H~ +t{wT,{u) + £) -u 

within the lattice ball of radius t^ with probability at least 1 — exp(— c^t). 
Here > as soon as e > 0. 

Corollary 1.1. For a p-perturbation tt of a locally regular CA and 
finite initial sets which fill space, 



in the Hausdorff metric, almost surely. 

Next is a generalization of the flat edge result [6] . In particular this implies 
that Lp — > Li when p — > 1, as promised. 

Proposition 1.2. Given a standard p-perturbation of a locally regular 
CA and any e > 0, there exists a p <1 close enough to 1 that Lp agrees with 
Li outside the e -neighborhood of the set of corners of Li. 

Our second theorem provides necessary and sufficient conditions for exact 
stability. Before its statement, it is instructive to look at the three super- 
critical Moore TG CA. The 9 = 1 case is additive and exact stability cannot 
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hold. (This can be proved by the methods of [6] or [25], but we give a differ- 
ent argument in Section 3.) For 6 = 2 one finds that Kij^ = co{M) and hence 
this is a quasi- additive case, that is, a CA with convex i^i/^. Quasi-additive 
CA share many properties with additive ones [10, 11, 13], and lack of exact 
stability turns out to be among them. Finally, in the = 3 case Ki/^ has 
16 vertices, of which three successive ones are (0,1), (1,2), (1,1), and the 
remaining 13 are then continued by symmetry. (This set, which the reader 
is invited to compute, is the innermost region of Figure 7.) Eight of these 
are the only points that Ki/^ shares with the boundary of its convex hull. 
In a sense, the fact that these eight vertices form a discrete set makes this 
CA as unlike a quasi-additive one as possible. This turns out to be precisely 
the condition needed for exact stability. 
Accordingly, we denote 

dK' = d{Ky^)nd{co{KyJ), 

and describe the relevance of properties of this set in our main result. 

Theorem 2. Consider a supercritical locally regular CA (which also 
satisfies our standing assumptions) given by T , with limiting shape Li, and 
its standard p-perturbation. There are three possibilities: 

Case 1. dK' consists of isolated points, no three of which are collinear. 
Then the following hold for p <1 close enough to 1; 

(51) Lp = Li. 

(52) Convergence to Li is tight: for any e > 0, there exists an M 
so that, for any t and x S tLi, P{x is within M of At) > 1 — e. 

(53) There exists a large C so that, with probability \, (t — Clogt)Lin 
7? C At eventually. 

Case 2. dK' consists of isolated points, three of which are collinear. 

Then (SI) and (S3) still hold for p <1 close enough to 1, but 
tightness (S2) no longer does. Instead, for any p <1 there exists a 
c > so that a corner oftLi is eventually at distance at least clogt 
from At, a.s. 
Case 3. dK' includes a line segment. 

Then (SI) no longer holds. Instead, for any p <1 there exists a 
c> so that a corner of tLi is at distance at least ct from At, a.s. 

Figure 1 shows a box neighborhood TG example for each of the three 
cases, from left to right, with periodic shading of updates: Case 1 (range 1, 
e = 3,p = 0.9), Case 2 (range 2,0 = 7,p = 0.95) and Case 3 (range 2,6 = 8, 
p = 0.95). 

The fundamental difference between Moore 6 = 2 and ^ = 3 TG CA is 
their mistake-fixing ability, which we now illustrate. Suppose we start each 
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case with a large copy of the invariant shape and remove a finite chunk of 
occupied sites at the boundary. Regardless of the location of such a hole, 
the 6 = 3 case eventually repairs (or "erodes") it and thus the hole's effect is 
bounded in time. Figure 2 provides a demonstration. This eroding property 
can be used to favorably compare the random dynamics on infinite wedges, 
determined by the corners of Li, to Toom rules [35]. The corners are then 
patched together by an oriented percolation comparison in the middles of 
the edges. In Case 2, mistakes are still fixed, but for wider wedges than in 
Case 1, and corners must be rounded off accordingly. 

By contrast, the = 2 TG CA can only repair holes away from the corners, 
while those at the corners have a lasting effect, as also seen in Figure 2. In 
a random dynamics, such mistakes pile up and induce a linear slowdown. 

Given the exact stability criterion of Theorem 2, it is natural to ask 
whether a typical supercritical CA has an exactly stable shape or not. As 




Fig. 1. The three cases of Theorem 2. 




Fig. 2. Error correcting for 9 = 2 and 9 = 3. 
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already mentioned, properties of typical monotone CA seem difficult to char- 
acterize. We will thus restrict our attention to a special family, TG CA with 
range p box neighborhoods Np. These are supercritical for < p{2p-\-l) [10]. 
The smallest examples are already illuminative. As A/i has already been dis- 
cussed, M2 is next in line and turns out to have = 1,2,3,5,8 in Case 3, 
9 = 7,9, 10 in Case 2, and = 4, 6 in Case 1. For very large ranges, ^'s in 
Case 3 form a small minority, as the following theorem demonstrates. 

Theorem 3. Fix an arbitrary e > 0. Among all supercritical range p box 
neighborhood TG CA, the proportion of those which are not exactly stable 
is for large p between 1/ log^'"''^ p and l/\og^ p. Here h = 2(1 — l/log2 — 
log log2/log2) f« 0.172. 

The proof of Theorem 3 connects the number of ^'s which lack exact sta- 
bility to the number of distinct products of pairs of natural numbers between 
1 and p. This latter is known as the Linnik-Vinogradov-Erdos problem, for 
which sharp asymptotic bounds were given by Hall and Tenenbaum [17]. We 
have no result on the division between Cases 1 and 2, but conjecture that 
Case 2 is much more prevalent. 

The rest of the paper is organized as follows. Section 2 contains the proof 
of a slightly weaker version of Theorem 1 and its Corollary 1.1. Section 3 
deals with Case 3, while Section 4 lays the geometric groundwork for the re- 
maining cases and proves Proposition 1.3. In Section 5 we introduce Toom's 
method and complete the proof of Theorem 2. Section 6 is devoted to a 
single example for which we can compute the shape for all values of the 
probability parameter p. In Section 7 we take a closer look at the collection 
of Kiiyj''s for fixed-range box neighborhoods, an analysis which culminates 
with the proof of Theorem 3. Finally, in Section 8 we finish the proof of 
Theorem 1 and discuss other related issues in lesser detail. 

2. Proof of Theorem 1. Recall that Theorem 1 deals with supercritical 
locally regular CA and their p-perturbations. These will be our context 
throughout this section. We will allow all constants C and c to depend on 
T and p in addition to their explicitly stated dependencies. (We emphasize 
that these constants will not, however, depend on the direction u.) In this 
section we only obtain a lower bound of the form 1 — exp(— Cef/log^t) on 
the probability of the event in Theorem 1. 

Many times below we will restart the random dynamics at a deterministic 
time or a random stopping time r. This simply means that only ^x,t with 
t'>T are used, with an initial state at time r which will be specified. 

Lemma 2.1. Assume that a finite Aq3{) fills the plane. Assume that x is 
at distance n from the origin. Then there exist constants c,C > ( depending 
on Aq) so that P{x ^ T^iAo)) < e'^'^ for k > Cn. 



10 



J. GRAVNER AND D. GRIFFEATH 



Proof. Call x surrounded at time t if x + Aq C Af. By supercritical- 
ity, there exists a time to at which ±ei and ±62 are all surrounded in the 
deterministic dynamics. Let Co = |T*o(j4o)|. If Po =P^°^ then itei and ±62 
are all surrounded at time to with probability at least po- Take a shortest 
lattice path p : = xq, xi, . . . ,Xn = x. We now define i.i.d. geometric(po) ran- 
dom variables Ti, . . . ,T„ as follows. Run the dynamics for time to- If a^i is 
surrounded at this time, Ti = 1, otherwise restart the dynamics with Aq at 
time to- Now run the restarted dynamics for time to; if it surrounds xi at 
this time, Ti = 2, otherwise restart again with Aq, and so on. In general, 
on the event {Ti = k}, Tj+i is the minimal i> 1 for which the dynamics 
restarted at time k + {i — l)to with Xi + Aq surrounds Xj+i at time to- 

By monotonicity and exponential Chebyshev, 

P{x i T'=*°(Ao)) < P{T^ + --- + Tn>k)< e-^'=^(exp(ATi))'^, 
for any A > 0. To conclude the proof, choose A small enough that -E(exp(A x 

ri))<oo. □ 

Note that this lemma implies that u'7r(n), if it indeed exists, is bounded 
away from uniformly in u, for any p > 0. 

Lemma 2.2. Assume that \T{Aq) \Aq\ = n, and start the p-perturbation 
from the same initial set Aq. // r = inf{t : T(Ao) C At}, then E{t) <p~^ x 
(log n + 3) . 

Proof. Note that all such sites attempt to get occupied simultaneously, 
each of them at each time with probability p. Hence r is geometric(p) for 
n = l. For n > 2, write a = — log(l — p) and divide the sum below into terms 
with k <a~^ logn and with k > a^^ logn to obtain 

00 

E{r) = J2{l-il-e-'^'r) 

k=0 

00 

< a'^ logn + 1 + ^(1 - (1 - e-"^n-^)") 

1=0 

00 

< a"^ log n + 1 - ^ n • log(l - e'^'^n'^) 

00 

<a-^logn + l + 2j2e~''' 
i=0 

= a-^ logn + 1 + 2(1 - e-")-^, 
and p=l — e^^ < a. □ 
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Proof of Theorem 1 with weaker probability estimate. With- 
out loss of generality, we can assume that u lies on or above y = \x\, that is, 
(n,e2) > l/\/2. Let = cr{^x,5 : s < t - 1, x € Z^}, t = 1,2,.. . . 

Let Tn be the first time (0, n) becomes occupied started from H~ and set 
Tn = Tn A Cn. By Lemma 2.1, -P(T„ ^ T„) < e"'^", for a large enough C and 
some c> 0. The crucial step is this bound: 

(2.1) \E{fn\Ts+i) - E{fn\Ts)\ < C'logn, 

for any s <Cn and some constant C . 

Recall that Tn is a deterministic function of (,x,t where {x,t) ranges over 
all space-time sites. As M is finite, Tn. depends only on a small subset of 
these variables. To be more precise, let comprise the sites {x^t) for which 
Tn depends on £^x,t- Then < Ctl^ and we can assume that the filtration 
ignores all other sites. At time s < Cn, let dAg consist of all the sites outside 
As which would become occupied if the deterministic dynamics were applied 
to As. Trivially, \dAs\ < |£„|. 

Restart the dynamics at time s + 1 with As. Let Tg be the waiting time 
after this at which all sites in dAs are occupied, that is, Ts = inf{k : OAs C 
^+i+fe}. By Lemma 2.2, E{ts\Ts) < C'logn. 

We now prove (2.1). We will repeatedly use the strong Markov property 
and monotonicity of the dynamics. To get the lower bound in (2.1), assume 
the worst case: no sites outside As (i.e., in dAs) get occupied, and therefore 
the dynamics faces an unchanged situation at time s + 1. Therefore, 

E{fn\J's+l)<EiTn\J's) + l- 

For the upper bound, assume that J^s+i reveals that all sites in OAs get 
occupied. Before we know J-s+i, we can only assume this happens after time 
Ts, and so the dynamics with the additional information is dominated by 
the one restarted at time s + Ts from the occupied set As U OAs- It follows 
that 

EifnlJ's) < E{fn\Ts+l) + E{Ts\Ts) < E{fn\Ts+l) + C'logn. 

This proves (2.1). 

Now let an = E{Tn), (in = E{Tn). By (2.1) and Azuma's inequality [19, 33], 

(2.2) P{\fn-an\ > s) < 2exp(-csV(^log^f^))- 
However, 

\an - «n| < £'(T„l|j'^>c„|), 

which is bounded by Lemma 2.1. From this it follows that 

P{\Tn - an\ >S)< P{\fn - a^j > s/2 -C)+ P{Tn -fn> s/2). 
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and after another application of Lemma 2.1 and suitable redefinition of c, 
(2.3) P{\Tn-an\ > s) < 2 exp{-cs'^/{nlog^n)) + 

For an integer i, let Ui be the largest j for which G . Then let 
be the first time at which all sites in B' = {{i, yi + n): \i\ < n^} are occupied. 
Moreover, let be the first time at which all the sites B" = : |i| < 

,yi + n — C < j < Ui + n} are occupied, where C is the diameter of the 
neighborhood M. Restart the dynamics at time with the occupied set 
at this time. Note that local regularity implies that within a constant time 
the deterministic dynamics occupies a large ball within a constant distance 
of any occupied point. By monotonicity, the deterministic dynamics would 
occupy B" in ti additional time steps, where ti is a constant which only 
depends on T. Applying Lemma 2.2 ti times one thus obtains 

EiT::-T:,)<C\ogn. 

Furthermore, let r^(^) be the time the dynamics reaches {i,n + yi) and let 
Tn{i) be the first time {i,n + yi) becomes occupied /rom the modified initial 
set 2/462 + H~ . The reason for this convoluted condition is that T„(i) with 
the same n are identically distributed, but this is not true for T^(i). 

To deal with different starting sets for T^(«), let Sn be the time the ran- 
dom dynamics fills H~ ni?(0,n^) from —62 + H~ (which is contained in all 
starting sets). By a similar argument as in the previous paragraph 

E{Sn)<C\ogn. 

Therefore, with a'^{i) = E{T^{i)), 

< a^(i) - a„ < E{Sn) < Clogn. 

Furthermore, the argument leading to (2.3) can be carried out with Tn re- 
placed by T^(i) and a„ by a'„(i). 
Therefore, for s > Clogn, 

P{Tn -T',>s)< P{\Tn - an\ > s/4) + J2 PiKii) - an{i)\ > s/4) 

|i|<n2 

< Cn^exp(-cs^/nlog^ n) + Cn^e~'''*. 

It follows that 

roc 

E{Tn - T'^) < CV^log^ n+ ^ P{Tn -T'^>s) ds, 

which after a short computation implies that E{Tn — T^) < C^/nlog^ n. 

We are almost done, but need an estimate for yet another approximation 
to Tn- Let Til' be the first occupation time of (0, n) started from B" — 1x62 = 
■ Nl <n^,yi-C <j < yi]. Then, for < < Cn, 

P{T';' -Tn>k)< P{T'^' + Tn) < P{Tn > Cn) < e-'^, 
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while for k>Cn, 

P{T'^' -Tn>k)< P{T;;' >k)< e"^'' 

by Lemma 2.1. Hence E{T^' — Tn) is bounded above by a constant. 

Now assume that < m < n. Restarting the growth process at time r„, 
we get 

am+n <a^ + an + E{T'^ - T„) + E{TZ - T^) < + a„ + C^^log^ n. 

By the deBruijn-Erdos subadditive theorem [33], a,„/n converges to a finite 
positive number a, which of course depends on p and u. We declare w-,^{u) = 
{u,e2)/a. 

To finish the proof, take first an outside H~ + ti(;^(ti)(l + e) -u. Let 
n = j-yi> twT,{u){l +e)/(n, 62) =t(l +e)/a. Then 

p((i,i)GAi) = p(r;(i)<t) 

<p(t;(.)<W(i + ^)) 

<Pi\TUi)-aUi)\>nae/2) 

< exp(— cn/log^n), 

for a large enough n. This proves the weaker version of the upper bound in 
Theorem 1. The lower bound is proved similarly. □ 

Several remarks are in order. First, note that the proof avoids the sub- 
additive ergodic theorem altogether, by combining properties of subadditive 
sequences with large deviation estimates. 

Second, it is in fact possible, by the same methods, to obtain a superad- 
ditive relation for an of the same order, namely, 

am+n > am + On - C^/nlog^ n. 

A closer look at the proof of the deBruijn-Erdos theorem (from [33]) then 
gives a rate of convergence for a„: |a„ — a| = 0(log^n/-^/n), which can be 
used to show that, within a lattice ball of radius t^. At is a.s. between 
{t±C^/t\og^t)■w■„{u)■u + H-. 

Third, the proof uses supercriticality and regularity only to "fill in." For 
any monotone, local, interface solidification with automatic coherence the 
proof remains valid. While we will not attempt to precisely define the con- 
cept, automatic coherence certainly holds when the interface moves upward 
(i.e., u = €2) and the growth is such that an empty site can never have an 
occupied site directly above it. Perhaps the simplest example is the ran- 
dom dynamics in which a site becomes occupied for sure with two or more 
occupied neighbors in its von Neumann neighborhood and with probabil- 
ity p with an occupied site directly below. Another class of examples are 
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the ET-exclusion processes [31]. For some of these examples, the fluctuation 
estimates mentioned above may be new. 

Finally, and curiously, there seems no way to make the proof work for 
general monotone dynamics which do not solidify. Such cases thus remain 
an intriguing challenge. 

Lemma 2.3. Fix an a > and e > 0. Then there exist constants c,C 
so that the following holds. Start the dynamics from Aq consisting of sites 
inside {H~ \ (—Cu + H~))riB{0,Cn). Then, An includes all sites inside 
B = {{H~ + nwT^{u){l — e)u) \ H~) n -6(0, Cn) with probability at least 1 — 

g— cn/ log'^ n 

Proof. Let T'{x) [resp. T{x)] be the first occupation time of x G B 
started from the stated (resp. from H~). By Lemma 2.1, P^sup^^^ T{x) > 
Cn) < e~'^". However, by a "speed of light" argument, on {supr^^^T{x) < 
Cn} the equality T(x) = T'{x) holds for all x £ B. The claim now follows 
from Theorem 1. □ 

Lemma 2.4. The function w-^ : 5*^ ^ R is continuous. 

Proof. Again assume that {u, 62) > l/\/2. For a fixed large C and small 
e > 0, H~ — Ctee2 C H~ C H~ + Ctee2, within the lattice ball of radius Ct, 
provided — f II2 < e/2. Let Ei,{k,t) be the event that all sites on the y-axis 
up to k are occupied at time t started from H~ . 

By Lemma 2.3 and Theorem 1, both the events £"^((1 — Ce)tt;^(u)t/(f , 62), t) 
and Ey{{l + e)wTt{v)t / {v .,€2) ,tY happen with probability (very) close to 1. 
This is only possible if (1 — Ce)w.,^{u) < (1 + e)wT^{v). An analogous reverse 
inequality is proved similarly. □ 

Proof of Corollary 1.1. An e > wih be fixed throughout this 
proof. 

For any direction u, Theorem 1 implies that with probability exponen- 
tially close to 1 



eventually. This is therefore true simultaneously for any finite collection of 
u's and then by Lemma 2.4, 



AtCH^ + tw^{u){l + e)-u. 



It follows that with probability 1 




— CH^ + Wn{u){l + e)-u, 



At 



C fl H- +w^{u){l+ef-u = {l + efK*y. 



t 
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eventually. 

For the lower bound, take a bounded, strictly convex, set D -^i/to^ (1 + 
e). Then = K* is and has for small enough 5 > the property described 
in the next paragraph. 

Start with consisting of sites inside nL^. Take k = jb Euclidean 
points xo,...,3;fe-i on the boundary of nL^, chosen so their directions are 
equidistant vectors in S*^, and let uq, . . . , Uk-\ be the outside normals to nLg 
at the chosen points. The enlarged set 

fc-i 

i=0 

includes {n + y/n)Li.. 

Now run the random dynamics from Aq for ^/n time steps. Since has 
boundary, we need to go just a constant distance inside to "see" the rele- 
vant portion of H~.. To be more precise, {{—Cui + H~.) \ {—2Cui + H~^)) n 
B{xi,Cy/n) is included in nL^, for all i. By Lemma 2.3, with probability at 
least 1 — exp(— C\/n/log^ n), all the sites in (n + ^Jn^L;, become occupied. 

Repeat the above procedure (running the random dynamics for ^Jn time 
steps) 3\/n times. As a result, {n^ j^/n)L^ C A^j^,-^ for j = 1, . . . , (in 

particular, 4nI/£ C ^4n), with probability at least pn = 1 — exp{—c^/n/ log^ n). 

Now fix an a < 1 and find a large fco so that nn>22''0 Pn> cl. Let Tq be 
the first time 2^'^"L£ C At. By what we proved so far, 

P((22fc + j2^-)L, c ^2fe+,-2^-To for i = 0, . . . , 3 • 2^ = 1, 2, . . . ) > a. 

We thus have a strictly increasing sequence of integers bm with bm+i — bm = 
o{bm), such that 

bmil-e)LsCAb,^ 

eventually, with probability at least a, thus a.s., as a was arbitrary. For any 
t between bm and bm+i, 

t{l - efLe C bm{l - e)Le C A„ C At 

eventually, finishing the proof of the lower bound. □ 

3. Lack of exact stability in Case 3. Fix a u G 5^. Let lu be the boundary 
line of —w{u) ■ u + . Note that w{u) is the largest number h>Q for which 
7r((— /i • u + ) n M) = 1 . Therefore, M Ciiu must contain at least one site. 

In general, for any line £ in the plane which does not go through the 
origin, let its open (resp. closed) lower cut C°{£) [resp. C~{£)] be the set of 
points in M which lie in the open (resp. closed) half-space of £^ which does 
not contain the origin. We emphasize here (as this convention will be used 



16 



J. GRAVNER AND D. GRIFFEATH 



extensively) that the points in C°{tj wih be called below the line i, and that 
left and right on the line are from the perspective of an observer who stands 
on i and looks toward the origin. 

We will make good use of duality between lines in Ki/^ and points of M in 
the sequel. The next lemma is our first example of this duality. To illustrate 
its statement (as well as the introduced terminology), let us consider an 
example. Assume that we are dealing with a TG CA and fix a direction u. 
Suppose also that iu contains an £ M such that a line i obtained by a 
small rotation of iu around has exactly 6 — 1 sites in C°{i). Note that for 
a sufficiently small rotation no other site but x^ is in £ n M. (An example 
for the range 2 Box TG CA with 6 = 8 is depicted on the right-hand side of 
Figure 3.) Therefore, for v close enough to u, is obtained by rotation of 
iu around Xu- A little geometric argument involving polar coordinates then 
shows that the boundary of i^i/^ must be flat at u/w{u). The lemma makes 
a stronger and more general statement and is illustrated by the left-hand 
side of Figure 3. 

Lemma 3.1. The following are equivalent for a . 

(1) There exists a line through u/w{u) which in a small neighborhood of 
u/w{u) lies in Ki/^;- 

(2) There exists a point Xu ^ iu^J^ so that if i is a line through Xu and is 
a rotation of iu by a small enough angle, TT{C°{i)) =0. 




Fig. 3. Illustration of Lemma 3.1 and its proof. 
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In case dKij^ is locally a line at u/w{u), Xu in (2) is unique. In fact, the 
smaller angle between dKij^ and u/w{u) is the same as the smaller angle 
between the vector and lu- 

Proof. Note that a short hne segment through u/w{u) perpendicular to 
the vector uq S is given in polar coordinates (with the angle represented 
by a unit vector v) by the collection of vectors 

{u,uo) V 

—-— ■- --.Wv-uW <a 

w[u) {v,uo) 

for a small a > 0. 

Assume first that the statement (2) holds. Let uq = —Xu/\\xu\\ be the unit 
vector pointing from Xu to the origin. Then (2) says that for v close enough 
to u, w{v) <w{u){v,uo)/{u,uq). It follows that 

(3.1) 4-> ^ 



w{v) w{u) {v,uo) 

The polar representation of a line mentioned above immediately demon- 
strates the implication (2) (1). 

To prove the reverse implication, note that (1) implies that (3.1) holds for 
some Uq, and let x'^ = —w{u)uo/{u,uq). Then x'^ has the properties required 
of Xu, except it may not lie in M. However, we can let Xu be the closest site 
in iu n AA. (We can in fact go in either direction from x^-) The fact that 
M is discrete ensures that a parallel translation from x'^ to x^ of any line i 
close to £u does not pass through any site of M. Thus (2) is satisfied. 

To prove the last statement, note that two different x^ would, by (3.1), 
produce two distinct open line segments, which would meet at u/w{u) and 
which would both be included i^i/^. But then a flat portion of dKij^ near 
u/w{u) would be impossible. □ 

Lemma 3.2. Fix a u £ which satisfies the condition of Lemma 3.1 
and pick a corresponding Xu- For v close enough to u, the concave wedge 
Q = II~ U H~ satisfies 'T{Q) C —Xu + Q- When dKij^ is locally a line at 
u/w{u), Q is invariant: 'T{Q) = —Xu + Q- 

Proof. The first part follows from Lemma 3.1: a point in T((5)n (— Xu + 
QY would imply that a point on the boundary of —x^ + Q sees a sufficient 
configuration in the interior of Q, but clearly —x^ is in the most advanta- 
geous position for this. This would translate, for £ as in (2) of Lemma 2.1, 
into 7r(£°(^) U £"(£«)) = 1, but if the rotation is sufficiently small (by dis- 
creteness of M) {C°{£) U C°{iu)) nAf = C°{1) nAf,a contradiction. 

The second part also follows because in this case 7r(£~(£)) = 1. For, other- 
wise Xu could be moved to the next point in ^„ n for which 7r(£~(^)) = 1. 
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[Again, such a point must exist or else w{u) could be decreased.] This would 
contradict uniqueness. □ 

When u satisfies the assumption of Lemma 3.2 there exists an invariant 
wedge of the following form: 

Q' = i-Mvi + H-)U {-Mv2 + H-) U H-, 

where vi and V2 are close to, but on different sides of, u and M is large 
enough. In particular, a hole of shape Q' dug into H~ may be translated 
by the dynamics, but is never filled. If the creation of such holes is random 
they pile up and, as we will demonstrate by the comparison process we now 
introduce, slow down the interface. 

The following randomly growing surface will be useful here and in Sec- 
tion 8. At every time t = 0, 1, 2, . . . a site x G Z has a height r/^(x) G Z+, 
with r/Q = 0. We will use two versions, which we call fast and slow, of the 
rule for increase in heights. Let b{x,t) be Bernoulli random variables with 
P{b{x,t) = 1) =p'. The slow version evolves according to the following rule: 

{ry^(x)-|-l, if 6(rE, t) = 1 and 
Vtiy) > Vtix) for ah y with |y - x| < 1, 
tjtix), otherwise, 

while the fast version updates as follows: 

{r/^(x)-|-l, if6(x,t) = lor 
Vtiy) > Vti^) foi' some y with \y — x\ < 1, 
r]t{x), otherwise. 

Note that the reverse dynamics, r]t{x) = t — r][[x) changes the version and 
replaces p' by 1 — p' . We will assume that b{x,t) are not necessarily in- 
dependent, but have finite range dependence in space: if either ti ^ t2 or 
|xi — X2I > r, then 6(xi,ti) and b{x2,t2) are independent. 

Lemma 3.3. 

(1) For the slow version: Given any p' > 0, there exist an a > and c > 
so that 

PiVtix) < at) < e^^*. 

(2) For the slow version: Given any e > 0, there exist a large enough p' and 
a c > so that 

P{v'^{x)<{l-e)t)<e-'\ 

(3) For the fast version: Given any e > 0, there exist a small enough p' and 
a c > so that 



P{r][{x) > et) < e" 
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Proof. The proof of (1) and (2) is a last passage percolation argu- 
ment. By [24] we can in fact assume that the random variables h{x,t) are 
independent. Once the neighborhood condition {'q[{y) > ri[{x) for all y with 
\y — x\< 1) is satisfied, a site x has to wait a geometric(p') number of time 
steps before it increases. Accordingly, let g{x, s) be i.i.d. geometric with suc- 
cess probability p' . By a simple inductive argument, it follows that the first 
time Tn{x) when r]'g{x) = n > 1 equals 

(n-l \ 
max< ^ g{xi,i) :Xn-i = x,Xi^i G {xj — l,Xi,Xi + 1} for < i < n — 1 >. 
I i=o J 

Hence 

Piv's (^) <n)= P{Tn{x) >s)< 3"P (j2 5(0' ^) > • 

By an elementary large deviation computation, we get, for a fixed p' > 
and a small enough a > 0, P{r]'g{x) < as) < exp(— cs), which implies (1). 
Another large deviation computation gives P{r]'^{x) < (1 — e)s) < exp(— cs) 
for a fixed e > and p' close enough to 1. Finally, (3) follows from (2) by 
reversal. □ 

Proof of Theorem 2 in Case 3. Let u £ be a direction of an 
interior point in a line segment of dK' . Then u/w{u) belongs to the interior 
of a line segment of dKij^ (satisfying the condition of Lemma 3.1) and the 
corner of L which corresponds to the edge of co{Kiiyj) containing u/w{u) 
moves, in the deterministic case, with speed w{u) in direction u. The follow- 
ing claim will therefore finish the proof. Start the random dynamics from 
Ao = i?t7 n Z"^. Then, for some a > 0, 

(3.2) P{At C (1 - a)w{u)tu + H~)>1- e'^^ 

For simplicity, rotate the space so that u = 62- Recall that 2M is the 
width of the bottom edge of Q' . Assign fn{i) = if Ai n (iMei +Q') = 0, 
and 7?i(«) = 1 otherwise. In general, let T]t{i) be the smallest k for which At H 
{iMei — kxu + Q') = 0- It is clear that fjt is for M large enough dominated 
hy r]t = t — rj[ where r][ is the slow version from Lemma 3.3 and p' = (1 —p)^, 
for some appropriately large k. The range of dependence r depends on M 
and angles between vi and u and V2 and u, but is clearly finite. Therefore, 
(3.2) follows from Lemma 3.3(1). □ 

4. Flat edges of shapes for p close to 1. The setup is the same as in the 
previous section. In Lemma 4.1 below, the direction of a rotation of a line 
i is determined by the direction of motion of the outward normal to the 
half-space in £^ which does not contain the origin. 
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The left-hand side of Figure 4 depicts a general situation in the statement 
and proof of the lemma, while the right-hand side again presents a TG CA 
example. This time the range 2 Box case has 6 = 1 (see Figure 6). Note that 
if ^ is a rotation of ^u around x^^ by a small negative angle, then there are 
exactly — 1 = 6 sites below I. The same is true for rotations around xj^ by 
a small positive angle. This translates to two line segments on the boundary 
of Kij^ which meet at u/w{u) at a convex angle, of 45° in this case. 

Lemma 4.1. For a u £ S^, assume that near u/w{u) the boundary of 
Ki/yj consists of two lines at interior angle below vr. Then there exist G 
iuCiM with the following properties. If i is a small rotation of iu either 
through x\^ by a negative angle, or through x^ by a positive angle, then 
TT{C~{i)) = 1. In fact, the smaller angles between dKi/^ and u/w{u) are 
the same as the smaller angles between and 1^ o.f^d between x^^ and if 
xf^ and xl^ are chosen to be furthest apart. 



Proof. The argument is very similar to the one for Lemma 3.1. The 
local equation for dKi/^ to the right of u/w{u) is given hy v ^ {v/uj{u)) ■ 
{{u,uo) /{v,uo)), for a suitably chosen uq. Then 




Fig. 4. Illustration of Lemma 4.1. 
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if V is to the left of u. If x'^ = —w{u)uq/{u,uq), this means that any small 
rotation £ of £u around x'j in the positive direction has 7r{C°{i)) = 1. Now 
simply move x'^ rightward on to the first point on ^„ H for which 
7r{C"{i)) = for small positive rotations i. Such a point must exist, or else 
7r(£°(^„)) = 1, which contradicts the definition of w{u). This defines xj^, 
which must be in MCil, as small rotations contain no other sites of M. The 
definition of is similar. □ 

Lemma 4.2. Fix u £ , and Lemma 4.1, chosen as far 

apart as possible. If vi is a small positive rotation of u, then the convex 
wedge = H~ n is invariant: 'T{W'^) = — xj^ + . Similarly, if V2 
is a small negative rotation of u, then the convex wedge W~ = n is 
invariant: T{W~) = — + W~ . 

Proof. This proof is completely analogous to the one for Lemma 3.2. 
We omit the details. □ 

Note that the two wedges from Lemma 4.2 are moving toward each other. 
In particular, if we now dig any finite hole in H~ , it gets filled, as we state 
more precisely in the next corollary. 

Corollary 4.3. Ifvi and V2 are as above, and 

Ao = {H- n {-Mv^ + H-)) U (H- n i-Mv2 + R-)), 

then 

T\Ao)=tw{u)u + H- 

for t > CM. 

Assume now that ui,U2 € are such that ui/w{ui) and U2/w{u2) are 
on the boundary of co^Kij^), and that U2 is a positive (counterclockwise) 
rotation of ui (by the smaller angle between them). Assume also that for 
small positive (resp. negative) rotations v of U2 (resp. ui), v/w{v) is not on 
the boundary of co{Kiiyj). This assumption always holds in Cases 2 and 3 
of Theorem 2 when ui/w{ui) and U2/w{u2) are vertices of co(-fCi/^), and 
also holds for other vectors in Case 2. Then 

{w{ui)ui + n {w{u2)u2 + R-^) = z + H~^n R-^ 

covers a vertex of L [and when ui/w{ui) and U2/w{u2) are vertices of 
co(/^i/^„) also portions of corresponding edges of L\. The equation above 
defines the vector z = z{ui,U2). While this wedge is by itself not necessarily 
invariant (although for small M it often is), a bounded perturbation with a 
suitably rounded corner is superinvariant [as in (3) of Lemma 4.5 below]. 
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Lemma 4.4. Let ui and U2 he as above. Apply Lemma 4.2 to u = ui to 
get the corresponding W~ and x^. The corner of W~ moves slower than 
IL^^, that is, {-xl_,U2) <w{u2). 

Proof. The conclusion is equivalent to x^^ ^ C~{£u2)- But this follows 
because xj^^ ^ C°{iu') for any u' between ui and U2- [If Ki/^ were a straight 
line between ui/w{ui) and U2/w{u2), xj^^ would belong to £u' for any such 
u'.] □ 

An analogous version of Lemma 4.4 of course also holds for u = U2 and 
the corresponding W~^. 

Lemma 4.5. Assume ui and U2 are as in Lemma 4.4. There exists a 
convex wedge Wui^u2 which is included in, and outside a bounded neighbor- 
hood of the corner equal to, H^_^ f] H^^ , such that the following properties 
hold. 

(1) When the part of the edge of co{Ki/^) between ui/w{ui) and U2/w{u2) 
is completely included in Kij^, 

(2) When the open part of the edge of co(i('x/^) between ui/w{ui) and 
U2/w{u2) has empty intersection with Ki/^, 

'^(Wnuu^) Z + {{Wu,,u2 + B2i0,a)) n Wu,,u2), 

for some a> 0. 

(3) In every other case, T{Wu-i,u2) ^ z + Wuj^^u2- 

Proof. The condition in (2) implies that every vector v strictly between 
ui and U2 has xj^^ G ^"{K) and x^^ ^ ^°{^v)- This, together with Lemma 4.2, 
readily proves (2): one simply makes Wu-^,U2 start with and end with wedges 
considered there, and connects them with a convex curve of small enough 
curvature. 

The condition of (1) implies that every vector v strictly between ui and 
U2 has x^^ G £y and x:^^ £ ^v- Again, Lemma 4.2 implies that there are a 
succession of invariant wedges which connect the rays with normals ui and 
U2- The final statement follows by a subdivision of the edge of co(-ftri/^^,) into 
subintervals of types considered in (1) or (2). □ 

Corollary 4.6. Fix an e > 0. There exists a convex set Lg C L, which 
agrees with L outside the e -neighborhood of its corners, so that for large 
enough M , 



T{M-Le) D {M + l)Ls. 
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Proof. The wedges from Lemma 4.5 can readily be combined to ap- 
proximate an arbitrarily large multiple of L, within an error confined to a 
constant distance from its corners. □ 

Arguably, oriented percolation is the most useful comparison model in 
random spatial processes. We now introduce the version we will use. While 
this is in fact a random perturbation of a one-dimensional CA, it is, as 
the name suggests, best to think about it as a random occupied set which 
tries to establish long-range connections. Sites (m,n) G Z_|_ x Z_|_ are either 
occupied or empty (n is often referred to as a level). The basic ingredients 
are Bernoulli(p') random variables b{m,n), m,n> 1, such that 6(mi,ni) 
and b{m2,n2) are independent whenever \mi — ni\ > r or ni 7^ n2. (It is 
important that r does not depend on p' .) Prescribe some occupied set in 
Z_|_ X {0}. For m > and n > 1, (m, n) is occupied if b{m, re) = 1 and at least 
one of its neighbors (m, n — 1) and (m — 1, re — 1) is occupied. 

Lemma 4.7. Fix any a G (0,1). Also fix a large integer M and let 
[0,M] X {0} be occupied. If p' is close enough to 1, then for large enough 
C = C{p') the probability of the following two events converges to 1 as 
M ^00. 

(1) Any {m,n) with an < m < (1 — a)n and re = 0, 1, . . . is within distance 
C log re of an occupied point. 

(2) For every re there is a connection (through neighbors) of occupied points 
from level re down to level which stays entirely in {{m,n) : re > 1, are < 
m < an + M + C log re} . 

Proof. These are standard applications of contour arguments (see, 
e.g., [5], so we omit the details. □ 

Proof of Proposition 1.2. l{yi,i = 0,. .. ,R—1, are vectors pointing 
to the R successive corners of L in a counterclockwise order, then 

(M + 1)L,= y (y^ + L,). 

We will now concentrate on the edge between yo and yi. Start the random 
dynamics from j4o = M • n Z^. Say that {m,n) G Z+ x Z+ is occupied if 
all the lattice sites in (re — m)yo + rreyi + ML^ are at time re included in An- 
If (m, re) is occupied, then both (m, re + 1) and (m -|- 1, re + 1) are occupied 
with probability at least p' which is a power of p given by the number 
of lattice points in (M -|- l)Le \ ML^. Moreover, given any configuration 
of occupied points on level re, points on level n + 1 which are r apart are 
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occupied independently. Here r is any integer such that (ryo + (M + 1)L^) n 
{ryi + iM + l)Le) = 0. 

The occupied points hence form an oriented percolation with a finite 
range of dependence. For p' close enough to 1, Lemma 4.7(1) shows that the 
conclusion holds with probability which converges to 1 as M ^ oo. However, 
this is enough as every set is covered in almost surely finite time. □ 

5. Exact stability in Cases 2 and 3. We continue with the setup of the 
last two sections. Let us begin with a statement of Toom's theorem. We 
will only state the version we need here (for which the proof is in [35]), 
although the conclusion holds in considerably greater generality [4, 36]. A 
two-dimensional Toom rule is a deterministic CA 7t given by a map ttt 
with the following property: 

(T) There exists a line £t which does not go through the origin such that 
TTTiS) = 1 if and only if either £°(^t) C 5 or A/" \ £-{1^) C S. 

Now introduce space-time error sites, those sites {x,t) for which b{x,t) = 
0. Here b{x, t) £ {0,1}, x £ , t = 1,2 . . . , are assigned before the dynamics 
starts. The state of the Toom rule with errors is then given by rjt E {0, 1}^ , 
which satisfies r]Q = l and 

7?t+i = TT{i]t) n {x : b{x, t + 1) = 1}, t = 0,l,.... 

To develop some intuition, note that without errors a finite island of O's in a 
sea of I's gets eroded by Ttj as it is "squeezed" between two half-spaces with 
boundaries parallel to £t- (However, this island may move in the process.) 
Thus (T) is often called the eroder condition. A natural question is what 
happens with such a rule under persistent introduction of low-density errors. 

Theorem 5.1. Ifr]t{x) = 0, then there exists a Toom graph G = G{x,t), 
whose vertex set is included in {{z, s) : s < t, z £ Z^} and which satisfies the 
following properties, for some sufficiently large C > 0: 

(1) The number of possible graphs G with m edges is at most C™. 

(2) For a graph with m edges there are at least m/G vertices which are error 
sites. 

(3) For any r > 0, if rit{y) = for ||x — y|| < r, then the number of edges of 
G is at least max{r/C, 3}. 

For the proof see [35] . 

In the classical application of Theorem 5.1, b{x,t) are i.i.d. Bernoulli(p). 
Then P{^t{x) = 1) converges to 1 as p— > 1, uniformly in {x,t). Thus £,t has 
an invariant measure with density close to 1. This also follows when b{x,t) 
are not independent, but have uniformly bounded range of dependence in 
spacetime. 
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Lemma 5.2. Assume ui andu2 are as in Lemma 4.5, with corresponding 
wedge W = Wux,u2 ■ Start a p-perturbation ofT from Aq = W nZ'^ . Consider 
the following two events: 

(1) Ex^t,M = is within distance M of At\. 

(2) F = {there is a C so that, within the lattice hall of radius t^ , {t — 
Clogt)z + W cAt}. 

Then for any e > 0, there exists an M so that for any point x £ tz + W the 
event Ex t m happens with probability l — ef uniformly in x and t ). Moreover, 

p{F)=i: 

Proof. It is convenient to translate the dynamics so that W is fixed, 
that is, consider A'-/. = (At — tz) n W. Also, rotate the lattice so that W has 
its maximum at the origin and ui and U2 are situated symmetrically with 
respect to the ?/-axis. It then follows from Corollary 4.3, Lemmas 4.4 and 
4.5 that there are finite constants C and to (which again only depend on T) 
so that the construction in the following paragraphs is possible. 

Cut a finite neighborhood of the origin with a horizontal line y = —C, and 
let to be the time the deterministic dynamics needs to fill W again if sites 
above the cut are removed. Run the random system in multiples of time to, 
with the proviso that if any site at time nt^ is above the cut, then all sites 
above the cut are set to immediately. Also make all the sites above the 
cut if during the time interval [(n — l)to + Ij'^io] a site within C of the 
sites of the cut does not become occupied because of a bad coin flip, that 
is, although the deterministic dynamics would make it occupied the random 
one does not. The resulting set of occupied points at time nt^ is called A'^. 

If an integer site x gW is not in A'^^, then either an integer site in W 
strictly below the horizontal line through x must be in ^n-i and a site on 
or above the horizontal line through x must also be in An-i, or else a site 
within distance C of x does not become occupied although the deterministic 
dynamics would make it occupied. In the latter case we call x an error site. 
It is clear that error sites have finite range of dependence in spacetime and 
occur with probability p' which is above a fixed power of p and thus can be 
made arbitrarily close to 1. 

By Theorem 5.1, uniformly for sites x G and n, 

P{x is at distance at least M from A'J < ^ C"'{p')'"/^ < {Cpf^'^. 

m>M/C 

(5.1) 

Now the claim concerning the event (1) readily follows. To prove that P{F) = 1, 
note that we can choose M = Clogn (this C does depend on p), so that the 
probability in (5.1) is below 1/n^ and thus 

P(some X gW r\ B{0,n'^) is at distance at least Clogn from A'^^) < C jr? . 
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From local regularity and Lemma 2.1, it now follows that 
P{{tz + W)n B{0,t^) <t At+cio^t) < c/t\ 
and Borel-Cantelli completes the proof. □ 

Proof of Theorem 2 in Cases 1 and 2. Construct a convex set 
in the following manner: 

Step 1. Start from nL, a multiple of the shape of the deterministic model. 

Step 2. Every corner of nL whose corresponding open edge of Kij^ contains 
directions u such that u/w{u) G dK' is "logarithmically rounded off" 
by introducing for each such u an edge with normal u and length 
Clogn (where C is some large constant whose value will become 
clear below), and ensuring that the resulting set is a convex subset 
of nL. 

Step 3. Round off every corner of the set obtained in Step 2 to produce 
locally a translate W = Wui,n2 of Lemma 4.5. 

If R is the number of directions u for which u/w{u) G (?(co(i^i/^)), then 
Step 3 has produced R wedges W , which we label Wq, . . . , Wr-i. Start with 
some large Lm and couple (by using the same £^x,t) the resulting R + 1 
dynamics: one started from integer sites in each wedge, and the last one 
started from those inside tL^j- 

If the percolation model introduced in the proof of Proposition 1.2 sur- 
vives for all time on each edge of Lt, then we call the coupling successful. 
This means that the state of each site is exactly the same as the state of the 
same site in one of the wedges. 

Fix an X G Lj. By the FKG inequality and Lemma 4.7, the following three 
events simultaneously happen with probability at least 1 — e if Af is large 
enough: 

(1) The coupling is successful. 

(2) Ex^t,M happens for a suitable wedge Wi guaranteed by (1). 

(3) F happens for every wedge Wi, i = 0, . . . , R — 1. 

For an arbitrary initial set which fills space, we once again use the fact 
that L]\f is covered in finite time to get (SI) (2) (3) in Case 1, and (SI) (3) 
in Case 2. 

It remains to show that the a.s. deviations from a corner in Case 2 are 
at least logarithmic. Pick a corner a £ L whose corresponding open edge of 
co(i^i/^^,) contains u/w{u), for some direction u. Note that the boundaries 
of nL and nw{u)u + H~ intersect at exactly ta. Finally, consider the infinite 
wedge W defined by the corner and locate its vertex at the origin. 

The number of sites in 



Sk = Wn i{-kw{u)u + H-)\{-ik- l)w{u)u + H-)) 
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is bounded above hy Ck. Let Tk be a geometric random variable with success 
probability qk = {I — p)^^''^ . 

Start the dynamics from sites in ML, where M is arbitrary. It is clear 
that 

P{ta is at distance at least ck from At) > P{Ti + • • • + < t) 

(5.2) 

(i) 

where Tf^ are i.i.d. copies of T^.. Now by Chebyshev 
(5.3) P(T,« + • • • + Tf > t) < '''"^^'^^-^ ^ 



it-kE{Tk)y qk{tqk-kY' 

and it is easy to check that when k = c log t for a sufficiently small c = c{p) 
the upper bound in (5.3) is 0{t~^/'^). The desired result now follows from 
(5.2) and Borel-CantelU. □ 

We note that logarithmic a.s. fluctuations (S3) are optimal in Case 1 of 
Theorem 2 as well, as any of the sites in (t + l)L\tL can stay unoccupied for 
time clogt as a result of bad coin flips. This happens independently for each 
such site with probability if c = c{p) is small enough. Since the number 
of such sites is linear in t, a large deviation computation shows that the 
probability that this happens for at least one site is at least 1 — exp {—cy/t ) 
[19]. 

6. An example. In this section we present an example of a one-parameter 
family of random rules TTp, p € [0, 1] , for which we can compute the half-space 
velocities explicitly. Every such example seems to be similarly based on the 
models introduced in [14] and [30]. Apart from the exactly stable cases, the 
example which follows therefore seems to be the only nontrivial instance of 
a random growth model with known shape. 

The best way to think about this model is on the hexagonal lattice, but 
we will describe it so that it fits into our setup. The model's neigh- 
borhood M consists of seven sites, the von Neumann neighborhood with 
two added diagonal sites: (1,-1) and (—1,1). Then vrp(S') = 1 when at 
least one of the following four conditions is satisfied: (—1,0) G S, (0, 1) G S, 
{(1,-1), (0,-1)} CS, {(-1,1), (1,0)} CS. On all other nonempty sets S, 
TTp{S) =p. This is a p-perturbation of the additive model, but not a stan- 
dard one. Nevertheless we will denote the half-space velocities by Wp and 
the shapes by Lp. Supercriticality and local regularity are trivial. 

Note that Tip interpolates between two supercritical growth models. When 
p = l, the CA is additive with neighborhood M, which thus has i^i/u,^ = M* 
[which is co(AA) rotated by 90°], and Li = co(AA). When p = 0, the vertices 
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of i^i/^o are (0,±1), (-1,1), (1,-1), (1,2) and (-1,-2), while Lq is a 
parallelogram with vertices at (±1,0), (—1/3,2/3) and (1/3,-2/3). What 
sets this model apart from a generic example is that certain initial sets make 
it exactly solvable for every p, in the sense that P{x € At) can be expressed 
as a Fredholm determinant of an explicitly known operator on l"^ [14]. These 
initial sets are four wedges, which together cover the plane: Wi = H^^ n HZ^^ , 
W2 = H- n H- , W3 = HZe, n H-^ , Wi = HZe, n HZ^, , where = (1, l)/\/2. 

Assume that Aq = Wi H Z'^. The first observation is that At = {{x,y) : 
y ^ gt{x)}, where gt'-Zi ^ [— oo,t] is a nondecreasing function. This is easily 
proved by induction. As a consequence, whenever x has its east (i.e., right) 
neighbor in At, it also has its southeast diagonal neighbor in At. On this ini- 
tial set, the rule is therefore symmetric across the line y = —x. Now let, for 
every positive integer n, ht{n) = gt{—t + n) — n. Then ht{n) = —00 for n < 0, 
ht{0) is a random walk which jumps by +1 (resp. stays put) with probability 
p (resp. 1 —p), and for n > ht+i{n) equals ht{n) + 1 = ht{n — 1) automati- 
cally whenever ht{n — 1) > ht{n). Finally, when ht{n — 1) < ht{n), /it_|_i(n) = 
ht{n) + 1 with probability p and otherwise ht+i{n) = ht{n) + 1. This estab- 
lishes the equivalence. It follows that there exists a self-invertible function 
: [p, 1] ^ [p, 1], so that At/t converges to the region = {{x, y) G : y < 
1, —(j){y) < X for p <y, and —l<x for y <p}. The function <f) has the fol- 
lowing explicit form [14]: 



The argument is similar when ^0 = ^2 H Z^. (In fact, that this case is equiv- 
alent to the above can be seen by mapping the model onto the hexagonal 
lattice, where it has fourfold symmetry.) If gt{x) = sup{— y + x £ At:y £Z} 
(with sup = —00), then ht{n) = gt{t — n) — n has the same evolution as the 
ht from the previous paragraph. It follows that this time At/t converges to 
Lw^ = {{x, y) € : y < 1, x < (j){y) — y for p <y, and x < 1 — y for y < p}. 
The remaining two wedge shapes are obtained by symmetry: = —Ly/^ 
and Ly/^ = —L\y2- 

The proof of the next proposition is very similar to the proof of Corol- 
lary 1.1, and hence omitted. (See also [10].) 

Proposition 6.1. Assume that a perturbation of a locally regular super- 
critical CA is given by vr. Assume that its initial set is a wedge: Aq = W , 
where W = D and ui and U2 form an angle in (0,7r). Then 



almost surely and in Hausdorff metric within any large ball of radius C . 




At 



t 
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This proposition, together with Corollary 1.1, immediately implies that 
in the present example, 



Hence the top half of the shape Lp is convex and C , but not strictly convex, 
for p> 1/2; is strictly convex and only at p = 1/2; and for p < 1/2 is 
strictly convex with a corner at its highest point (— yo/2,yo) above the x- 
axis. See Figure 5 for a plot. 

We note that fluctuations from the limiting shape in every direction, ex- 
cept ±(—1/2,1) when p< 1/2, can be obtained from [14]. For example, 
consider a e (—1, {—yo/2) A {—p)) and let gt{x) = in{{y eZ:{x,y)e At}. (It 
is easy to see that sites in At with a fixed x-coordinate always form an in- 
terval if they do so at t = 0.) Then (5t([aiJ) — </)(— a)t)/t^/^ converges in 
distribution to a nondegenerate random variable. This follows from the fact 
that for such a the evolution of At starting from a finite set and from Wi H 
can be coupled so that the difference of respective gt{lat\) is stochastically 
bounded. 

7. Exact stability for box neighborhood TG CA. The TG CA with box 
neighborhood of radius p has p{2p+ 1) supercritical thresholds 9 [10], and the 
same number of corresponding Ki/^ which we label Ki, . . . , Kp(^2p+i), and 

superimposed with different shades in Figure 6. Let £ = £p = U9=r^^ dKg. 
At first this set appears to be of bewildering complexity (cf. the range 5 
example on the right-hand side of Figure 6). However, perhaps the first 
feature revealed upon closer inspection is that £ consists entirely of straight 
lines, called K -lines, which extend through the entire picture. [In fact, if we 
were to include the critical Ki/^ for 9 = p{2p + 1) + 1, . . . , p{2p + 1) + p all 
these lines would continue indefinitely.] There is one-to-one correspondence 
between the ET-lines and the points of AA\ {0}, as we now explain. 

For any of the {2p + 1)^ — 1 sites x £ J\f\ {0}, start with a line through 
X and 0. Then rotate this line in the positive direction until it hits again. 
Call all such rotations ix,(j>, < (p < 2tt. The set of all cardinalities A(j;, (j)) = 
\C^{lx,<i))\ is exactly the set of 9 for which x S for some u. Moreover, 



and therefore its top half comprises points (x,y) which satisfy 

-(t){y) <x<(t){y) - y, i(p<y<yo, 
— l<x<l — y, ifO<y<p, 



where 




P>l/2, 
p< 1/2. 
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Fig. 5. Shapes Lp for p = 0, 0.1, ... ,1. 

by Lemma 3.1, whenever x is the only site in ix,(j)C^-^^ this hne determines 
a direction pointing toward the interior of an edge of K^(^^^^^, namely the 
direction of the normal to ix^cf, which points toward the origin. When the 
number of points in ix,<j} n to the left and to the right of x are equal, the 
direction points to the interior of an edge in Kg for 9 = A{x,(j)) — | {points 
to the left of x}| . 

Similarly, a line containing exactly two sites in AA \ {0} determines a 
direction in which exactly two dKg meet (at a point which is a vertex of 
both, for one a convex vertex, for the other a concave one). A line containing 
exactly three points in AA\ {0} determines a direction in which exactly three 
dKo meet (but this point is a vertex of only two of them). And so on. We 
are now ready to prove the next proposition, which in particular allows 
unambiguous reconstruction of all Kg from £. 

Proposition 7.1. Two different dKg intersect only in a discrete set of 
points. Moreover, all finite tiles of £ are triangles or quadrilaterals. 

Proof. The first statement follows since for all but a discrete set of 
rotations (j), x is the only site in lx,<j> H M and so the corresponding edge 
lies only in Kf^(^x,<f>)- Fix a ^-nd assume ^x,(j) = for some direction u. If 
£u contains two or more points in M, then dKg^ does not intersect dKg^+i 
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Fig. 6. The 10 (resp. 55) supercritical Ki/^,'s for range 2 (resp. range 5). 



if and only if \C~{£u)\ = Oo- Let x be the rightmost point in n such 
that ix,cf> = foi^ some (p. Decrease cf) to the largest cj)' such that £x,4>' H M 
contains at least two points. Clearly £x,(j>' n is a line for ^ = + 1 (and 
perhaps some larger 0's). A similar argument holds in the other direction, 
proving that, among two consecutive vertices of dKg^, at least one belongs 
to OKq^j^i. This is clearly enough. □ 

For any x, let A*{x) =mi^\A(x,(j))\- 

Proposition 7.2. dK' includes a line segment if and only if 9 = A*{x), 
for some x. Moreover, for every x, A*(x) = |A(x,0)| for some (j) such that 
£x,tj} exits Af = Boo(0,p) C through the neighboring (rather than opposite) 
sides. Finally, each K-line includes exactly one line segment on dK' , for 
some 9. 

Proof. Pick an x and (j)Q so that x is the only site in ix,4,Q H AA, and 
let 9 = A{x,4>q). For this 9 and normal u to £x,<f>Q-, = £x,<f>o- For the first 
assertion, it suffices to show that \i 9 > A*{x), then u/w{u) cannot lie on 
d{co{Kg)). But if it would lie there, all of i^i/^ would have to lie on one 
side of the line of Ki/^ determined by x and (po (see Lemma 3.1), meaning 
that A{x, (j)) >9 for every (p, a contradiction. 

For the second assertion, assume the said line exits the left- and right- 
hand sides of the square. We can also assume that x lies either strictly inside 
the third quadrant, or on the negative y-axis. In both cases, rotate the line 
around x to angle (p' in the negative direction just past the southeast corner 
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of the square; this results in A{x,(j)') < A{x,(j)) (with equahty in the second 
case) . 

For the final claim, assume that x is as in the above paragraph and that 
the line £x,(i} produces A*(x). Now rotate it in the negative direction to the 
smallest angle (j)' for which the number of points on the left of x on ix^c/,' CiM 
is larger than the number of points on the right. If a and b are the lengths of 
the line segments on ix^tj,' from the left and bottom edges of M, respectively, 
then a> b. This remains the case for any <j)" < (p' and thus further rotations 
in the negative direction only lose more points. An analogous argument 
works for positive rotations. □ 

Therefore, each K-line contributes exactly one edge on exactly one dK' . 
What produces a prevalence of exactly stable cases are those dK' which have 
more than their share of edges. For example, it is immediate by symmetry 
that the number of edges on any dK' is either or at least 4. The number of 
9 with lack of exact stability hence does not exceed ((2p+ 1)^ — 1)/4 = + p, 
and therefore the number of exactly stable ones is at least p'^ . This argument 
is quite simple, yet it fails to produce any way to identify a single exactly 
stable case. Our next result remedies this somewhat. However, we do not 
have an algorithm which lists more than 0{p) cases of either type. 

Proposition 7.3. All 6>2p^ + 1 are exactly stable cases. On the other 
hand, all 9 < p and 6 = p + 1 + i{2p +1), i = 0, . . . , p — 1, are not exactly 
stable. 

Proof. For the second assertion, consider first 9 < p and take x = (— />+ 
9 — 1, —p). For any angle 0, A(x, (j)) > 9. This immediately implies that the 
edges in Ki/^ adjacent to the positive y-axis lie in dK' . If 9 = p+1 + i{2p + 
1), i = 0, . . . , p, there is instead a single edge, perpendicular to the y-axis. 

For the first assertion, it is enough to prove that all first-quadrant bound- 
ary edges of such Kg have slopes in [1, oo]. Equivalently, take x £ M in the 
third quadrant strictly above the line through the origin with slope 1, and 
a line i through x with normal a negative rotation of 62 by angle at most 
7r/4 and such that x is the only point in iCiM. Then < 9. This is 

certainly true when u is close to vertical, and further rotations can only 
decrease □ 

Note that the proof of the above proposition shows that the 9 = 2, . . . , p 
thresholds have at least eight edges in dK' , which improves the lower bound 
on the number of exactly stable cases to p^ + p-1. This is a good lower 
bound for small p, although the exact enumerations are taxing. 

Proof of Theorem 3. By Proposition 7.2, we can reformulate the 
problem as follows. Consider all integer points {a,b), l<a,b< p. Take a 
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line i through (a, b) which intersects the positive halves of both axes. For 
each such line, let A(^) be the number of integer points in the closed triangle 
T{£) C bounded by £ and the positive halves of the axes. Finally, let 
i* = i*(a, b) be a line which minimizes T{i). We need to find an upper bound 
for the number of different T{i) over all (a, b). In the following computations, 
0(1) refers to a term of arbitrary sign whose absolute value can be bounded 
by a constant independent of p. 

Let i* and i* be the line segments on i* from (a, 6) to the x- and y- 
axes, respectively. The first observation is that £* can be chosen so that 
the lengths of and i* differ by 0{1). (In fact, these lengths can be made 
arbitrarily close, although not necessarily equal.) In particular, the area of 
T{£) is 2ab + 0{l). We will assume, without loss of generality, that the length 
of i* is not larger than the length of . 

Now find an integer point (0, yo) on the y-axis immediately below where 
£* intersects the y-axis. Reflect (0,yo) through (a, 6) to get a point {xi,yi) 
within 0(1) of the intersection of ^* with the x-axis. [Note that (xi,yi) lies 
within the closed first quadrant.] Also, let (xq, 0) be an integer point on the x- 
axis immediately to the left of where i* intersects it. Form the closed polygon 
n C by connecting (a, 6) (0,yo) ^ (0,0) (xo,0) {xi,yi) ^ {a,b). 

The area of 11 is 2a5 + 0(l). Most importantly, the number of integer 
sites in IT is A(£*) + n*/2 + 0(1), where n* is the number of integers in dli 
which are not on the axes. This follows since one loses about half of these by 
a small rotation of the line from (0,yo) to (xi,yi) around (a, 6). Therefore, 

n* 

T{t) = |{integer points in the interior of n}| + y + 2a + 26 + 0(1), 
and by Pick's theorem [1], 

n* 

area of 11 = | {integer points in the interior of 11} | + — + a + 6 + 0(l) 

= T[t)-a-b + 0{l). 

It follows that 

T{t) = 2ab + a + b + 0(1) = \{2a + 1)(26 + 1) + 0(1). 

Let Mat be the number of different products mn of integers m,nG [1, N]. It 
follows that 

|{r(r ) ■.o<a,b<p}\<c- M(2p+i) < c • 

log"' p 

by the Hall-Tenenbaum sharpening of a theorem of Erdos ([17], Theorem 
23). This ends the proof of the upper bound. The lower bound follows be- 
cause the lower bound in [17], Theorem 23 is obtained using odd integers 
only. □ 
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8. Final remarks. In this section we mention several results which are 
related to the main topics of the paper, sketch their proofs, and also complete 
the proof of Theorem 1 by removing log^ t from the large deviation estimate. 

Remark 1 (Continuous time). A standard continuous-time growth model 
At is obtained by adjoining every site x at an independent rate 1 exponential 
time after the time Tx at which x £ T{At^). This process can be constructed 
in the standard way by attaching a Poisson process to every x. Theo- 
rem 1 is still valid in this case. The a priori large deviation bound however 
has log^t replaced by log^i. We now sketch the proof. 

Observe At in discrete time units t= 1,2, Change it to A[ by making 

sure that between each time t and t + 1 no site at distance more than C log t 
from A't gets occupied. As is easy to see by comparison with the continuous 
time additive dynamics having the same neighborhood, 

P{At 7^ A't within a lattice ball of radius t^) < t~^, 

when C is large enough. Now continue the proof with A't, which of course is 
a discrete-time monotone Markov process. Lemma 2.2 must be used Clogt 
successive times to obtain the analog of estimate (2.2), and here is where 
the larger power of log originates. From this point the proof proceeds on 
familiar grounds, yielding existence of the asymptotic speeds w, while the 
continuous-time version of Corollary 1.1 establishes existence of the shape L. 

Remark 2 (Approximating half-space velocities) . Perhaps the most con- 
venient method for simulating a growth CA started from a half-space is to 
use a strip with tilted periodic boundary. Take a vector u at angle (f) £ [0, 7r/4] 
to 62. (By rotations and reflections it is clearly enough to consider these.) 
Then take a large L, and restrict the growth to the strip Hm = [0, M — 1] x Z. 
Let K = tan(f). Given any configuration of occupied sites inside Hm, extend 
it to Hm = [—M,2M — 1] x Z, by identifying the state of (x,y) with that of 
(x - M, [{y - kM)\ ) if X > M and with that of (x + M, [{y + kM)\ ) if x < 0. 

Start from Aq = {{x, y) G Hm ■ y < kx}. Let the dynamics update sites in 
Hm with the specified boundary conditions until some large time t when the 
interface apparently equilibrates. At this point, the average height above all 
points in [0, M — 1], multiplied by cos(j)/t, is a good approximation to Wt^{u). 

Note that this is a much more efficient technique for computing the shape 
Lp than merely running the dynamics from a finite seed and observing 
the resulting blob. In particular, smoothness of Lp is impossible to discern 
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that way. The method outUned above, on the other hand, uses averaging 
to greatly reduce transversal fluctuations on the interface. The theoretical 
underpinning is partly given in our last theorem. 
For a fixed M, and any x £ [0, M — 1], let 

hl{x) = max{y : {x, [y — kx\ ) G At}, h\ = max{/it(x) : a; G [0, M — 1]}, 

=min{y : (x, [y - kxJ) ^ At}, hP^ = m\n{ht{x) : x £ [0,M - 1]}. 



Theorem 4. Fix an e >0. If M is large enough, then with probability 1 
— < lim mi — < iim sup — < — , 

COSCp i->oo t t^oo t COSCp 

uniformly in u. 

In fact, it is easy to show by subadditivity that ht/t and hj/t both con- 
verge a.s. to the same number as t — > oo. 

Proof of Theorem 4. We start by proving the lower bound. Note that 
the boundary eff'ects spread with finite speed, as M is finite. Thus, until the 
time to = cM, the occupied sites on any vertical line through x S [0, M — 1] 
are above those started from an infinite tilted half-plane through (x, [kx\) 
or through (x, [kx\ — 1). (We have to allow for the second possibility because 
there may not be a perfect match at the boundaries.) By the weaker form of 
Theorem 1, the probability that the lowest unoccupied site above a fixed x 
is below Kx + {w.,r{u){l — e)/ cos (j))to is at most exp(— cto/log^ to). For n > 0, 
define this translation of Ao'- 

„ f / \ Wt,-(u) — e } 

Bn = <{x,y):y < KX -\ — nto > 

[ cos (p ) 

and the event 

E = {BiCAt,}. 

It follows that 

P^E") < Mexp(-cto/log^to) < e, 

if M is large enough. 

Now run the dynamics until time to- If E happens, restart the dynamics 
from the set Bi; otherwise restart the dynamics from Bo = Aq- Then repeat 
from the possibly translated Aq. Let Un be the largest u for which Bu C Anto ■ 
We have just proved that f/„ dominates an n-step random walk which at 
each step increases by 1 with probability 1 — e and stays put with probability 
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£. Therefore ?7„ > (1 — 2e)n with probability at least 1 — exp(— cn). This ends 
the proof of the lower bound (as uniformity in u follows because the constant 
Ce in the weaker form of Theorem 1 does not depend on u). 

The upper bound is proved similarly, except for the fact that we need 
an upper bound on the extent to which At can propagate in Iq steps. The 
trivial bound Cto, where C is the diameter of N , suffices. The comparison 
random walk now increases by 1 with probability at least 1 — e and by C 
with probability e, and so its speed is 1 + Ce. □ 



Theorem 4 remains valid for continuous-time growth as well, with a similar 
proof, except for two significant differences. The first is that the coupling 
of finite and infinite systems only holds up to time cL with probability 
exponentially close to 1 in L [15]. The second is that the trivial upper bound 
at end of the proof is not available, so the forward jump of the comparison 
random walk is arbitrarily long, with exponential tail probabilities. It is clear 
that the proof is still valid under these conditions. 

Proof of Theorem 1 (concluded). We choose u, M and io as in the 
proof of Theorem 4. Recall that to/M is small and so during a time interval 
of length to sites at a distance larger than M cannot interact. Define 

( , s , ■. WTr(u) — £ 1 

Bni = { (x, y) : iM <x<{i + l)M - 1 and y < + ^ ' nio \ 

L coscp J 

and the event 

E = {Bi^Q(ZAt,}. 

As before, P{E'^) < e \i M is large enough. Now we define rjuii) to be the 
maximal k such that every site {x,y) with x G [iM, [i + 1)M — 1] and y < 

"'cos!/) ^ ^^0 + '^^ ™ ^nto- We can couple Ty„ and the slow version of rj'^ 
from Lemma 3.3 (with p' = 1 — e), so that fjn{i) > ??^i(i) for every i and n. 
(The Bernoulli random variables b are probabilities of suitable translates of 
the event E.) It follows that, for every 6 > we can find a small enough e 
(which then dictates a large enough M), so that 

P(i0,y) G At for all y < (1 - > i _ g-c*. 

V COS (p ) 

This proves the exponential bound on the probability that Atjt lags signifi- 
cantly behind w-,^ (u) . We will now show that it cannot progress significantly 
faster than Wt^{u) either. To this end, we redefine 



Bn^i = ^{x,y) -AM < X < {i + l)M - 1 and y < kx + 



— nio 

cos 
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and 

E = {At,ClBifl}, 

so that again P{E^) < e. In case E fails, the occupied sites above the interval 
[0, M — 1] cannot progress by more than the diameter of M. Furthermore, 
now 'fjn{i) is the maximal k such that some site {x,y) with x G [iM,{i + 
1)M — 1] and y < Rk + ^ cis 0'"'^ + ™ ^"*o • Here is a suitably large 
multiple of the diameter of M, which ensures that the fast version of r]'^ in 
Lemma 3.3 (now with p' = e) dominates fjn- Thus Lemma 3.3(3) completes 
the proof. □ 

Remark 3 (Continuity of Ki/^^). Assume a standard p-perturbation of 
T. As p changes from to 1, Ki/wp varies continuously. To see this, assume 
that p' is close to p, p' <p and couple the systems with the two probabilities 
in the obvious way. To show that Wp'{u) is close to Wp{u) (uniformly in ti), 
one needs to look at the proof of the lower bound in Theorem 4. Between 
and tlie occupied sets in the two systems will not differ at all with 
probability (1 — {p — p'))^*" ■, which, as to is constant (albeit dependent on e), 
can be made larger than 1 — e \ip — p' is small enough. Once this observation 
is made, it is only necessary to follow the rest of the lower bound proof with 
p replaced by p' . 

Note that this continuity alone demonstrates that Lp has corners (i.e., is 
not differentiable) in every case not quasi-additive when p is close enough 
to 1 (although these corners may not move at the same speed, or even in 
the same direction, as the corresponding corners of Li). 

Remark 4 (Shapes for small p). Again, assume a standard p-perturbation 
of T. What happens as p ^ 0? Certainly Lp shrinks, and in fact 

—Lp > L, 
p 

the limit shape of the continuous-time growth model At. To see this, let 
A'^ be Ayt/pl- After a site sees a sufficient configuration in A[ (resp. At)., it 
becomes occupied in a time distributed as Tg (resp. Tg). It is easy to see 
that for small p, distributionally, Tg > T(,{'i — p)- Rescaling of continuous 
time immediately gives L/[l — p) D Lp/p for small p. 

For the opposite direction, observe that Te+p>Tg, in distribution. The 
lower bound part of the proof of Theorem 4 for continuous time now shows 
that AtQ D B\ with probability 1 — With the same probability, then, A'^ D 
Bi at time t = to + CpLto < (1 + s)to if p is small enough. (The added term 
is simply p times the number of sites in Bi \ Aq.) This easily completes the 
proof. 
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In closing, let us pose two challenging conjectures based on experiment. 

Conjecture 8.1. If tt is a standard p-perturbation of a locally regular 
and supercritical CA with convex Ki/^, then Kij^^ is strictly convex. 

This conjecture fails for nonstandard perturbations, as seen from the ex- 
ample discussed in Section 6. In fact, it fails for nonstandard perturbations 
even if we restrict to p very close to 1. A range 2 box counterexample is 
obtained by '/r(5) which is 1 when |5| > 9 and p when |5| = 8. A glance at 
Figure 6, together with results from Section 3, confirms that Kijyj^ cannot 
be convex for any p <1. 

Figure 7 illustrates the application of Theorem 4 to two examples. The left 
frame depicts Kij^^ for the Moore TG CA with 9 = 2> and p= 1, 0.9, . . . , 0.4, 
while the right frame does the same for the range 2 box TC CA with 
9 = 9 and p = 1, 0.975, . . . , 0.5. As guaranteed by Theorem 2, co{Ki/^^) C 
co(-ftri/^^,^) for p close enough to 1. What is more, the angles at the cor- 
ners of Ki/wp approach the angles of Ki/^-^ as p^l. (This can in fact be 
proved by methods of the present paper.) As p decreases, dKi/^^ separates 
from dKijyj^, but Ki/wj, remains nonconvex. (Whether it may become 
before the boundaries separate is not clear.) Upon further decrease in p one 
observes concavities gradually filling in, until Ki^^^ becomes convex. Such 
observations, as well as early belief in the asymptotic isotropy of Eden's 
continuous-time random growth model [7], suggest our final conjecture. 

Conjecture 8.2. If p is small enough, the standard p-perturbation of 
T has strictly convex and smooth K^^^^. The continuous-time version K\i^ 
is also strictly convex and smooth. 




Fig. 7. Two examples of TG Ki/^^'s. 
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